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ABSTRACT 

The contribution of the quiet Sun to solar irradiance variabihty, due either to 
changes of the solar effective temperature or to the presence of unresolved mag- 
netic field, is still poorly understood. In this study we investigate spectral line 
diagnostics that are sensitive to both temperature variations and the presence of 
small scale unresolved magnetic features in these areas of the solar atmosphere. 
Specifically we study the dependence on the magnetic flux density of three pa- 
rameters describing the shape of two magnetically sensitive Fe I hues, at 630.2 
nm and 617.3 nm, namely the line core intensity (IC), full width at half maximum 
(FWHM), and the equivalent width (EQW). To this aim we analyze observations 
of active region NOAA 11172, acquired with IBIS at the Dunn Solar Telescope, 
as well as results from numerical synthesis. Our results show that IC is sensitive 
to both temperature and magnetic flux density variations, FWHM is mostly af- 
fected by magnetic field changes, and EQW is mostly sensitive to temperature. 
Variations of a few percent of the measured line parameters are found in ob- 
servational data that was spatially degraded to represent quiet-Sun, disk-centre, 
medium resolution observations. It is therefore possible to disentangle magnetic 
from pure thermodynamic effects by comparison of temporal variations of the 
EQW and the FWHM of either the two Fe I hues. 

Subject headings: Sun: photosphere - Sun: surface magnetism - Radiative transfer 
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Introduction 



It is well established that both the total and spectral Solar Irradiance (SI, henceforth) 
are modulated by the emergence and evolution of the magnetic field in the solar atmosphere. 
In particular, SI variations measured on time scales from days to several solar rotations 
can be explained by accounting for the radiati ve effect of the presen ce of magnetic flux 



in sunspots, faculae, and network regions (e.g. 



Domingo et al. 



20091 ) . However, there is 



still uncert ainty in measurements of SI variat ions on time scales of the activity cycle and 



longer (e.g. 



Harder et al 



2009 



Frohlich 



201ll ). as well as in the determinati on of individua l 



contr i butions of different m agnetic features to the measured variations (e.g. 



2011 



Krivova et al. 



Fontenla et al. 



20121 ). The role of the quiet Sun, that is those areas on the solar 
surface away from magnetically active and network regions, the so called Inter Network 
regions, is particularly debated. For instance, long-term SI variations unexplained by some 
recent models have been at tributed either to a 0.2 K global change in the Sun's effective 



temperature ((Frohlich 



20111 ) or to quiet Sun variations assoc iated with small-scale magnetic 



field changes during the activity cycle (IFontenla et al. 



2OI2I). 



Two main difficulties hinder the verification of these hypotheses. First, investigations of 
temperature variations would be achievable only by measurements of absolute radiometric 
contribution of the quiet Sun, and such measurements are not available yet. Second, the 
small-size magnetic features (10 km) that are se en permeating the whole solar surface 



in high-resolution observations (IBonet et al. 



2012 



and references therein), are unresolv ed 



in full-disk observations usually employed for SI reconstructions ( iDomingo et al. 



20091 ) 



so that their contribution to SI variations is hidden. Note that the origin of this so 
called quiet-Sun magnetic field is also controversial, as it is not clear yet whether and 
to what e xtent, it re s ults from the decay of active regio ns, turbulent dynamo effects or 



both (e.g. 



Lites 



2011 



Lopez Ariste fc Sainz Dalda 



2OI2I ). Variations of quiet Sun magnetic 
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flux over the cycle would indicate a global dynamo origin of the field; nevertheless, 
most of such measurements seem to indicate no temporal variations of quiet Sun field 

, and references therein) thus suggesting a 



( ISanchez Almeida Sz Martmez Gonzalez 



2011 



local dynamo origin. Whether this result is due to lack of sensitivity of instrumentation and 
the observational techniques employed, or to a true constancy of the field over the cycle is 
still an open issue. 

The presence of magnetic field influences the emitted radiation in magnetically sensitive 
lines directly, by the Zeeman effect, and indirect ly, by changing the thermodynamic 
properties of the plasma (e.g. 



Fabbian et al 



20101). Variations of line properties have been 



there 



Stenflo fc Lindegren 



'ore employed to study the properties of quiet Sun. For instance 
(1l977l ) derived, by the investigation of the subtle broadening of a large sample of Fe I lines, 
an upper limit for the magnetic field amplitude of quiet Sun network of 1000-2000 Gauss 
and an upper limit of 90 Gauss for intranetwork regio ns. Moreover, liiie para meters of a set 



of Fraunhofer lines were observed from 1974 to 2006 (ILivingston et al 



20071 ). with the aim 



of inferring long-term variations of th e photospheric temperature. Analysis of these data 



lead, however, to inconclusive results ( Penza et al. 



2004 



Livingston et al. 



2003). 



In this study we investigate the effect of magnetic features on the shape of two 
Fe I magnetically sensitive photospheric lines, the 630.2 and 617.3 nm, by analysing 
high-resolution spectro-polarimetric data of NOAA 11172 observed on March 17th, 2011 
at the Dunn Solar Telescope with the Inte rferometric Bidimensional Spectrometer (IBIS 



Cavallini 



2006 



Reardon &: Cavallin 



20081 ). Qualitative comparison with results from 
numerical synthesis has been performed to disentangle between the direct and indirect 
effects of magnetic field on the measured line shapes variations. We then investigate the 
effects of unresolved magnetic features on the shapes of the two lines by analysing IBIS data 
spatially degraded to mimic medium resolution observations. The Fe I 630.2 and 617.3 nm 
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lines were chosen because they have been frequently employed for photospheric magnetic 
field studies. Moreover, data at these spectral ranges are currently acquir ed with sever al 



instruments, like the E 



2003|), and SDO/HMI (IScherrer et al. 



INODE/SOT 



Tsuneta et al. 



2012). 



2008|), SOLIS/VSM flKeller et al. 



The article is organized as follows: in Sec. 2 we describe the observations, the data 
reduction and the numerical synthesis; in Sec. 3 we discuss the line formation properties of 
the two lines; in Sec. 4 we present the results derived from analyses of observational data; 
in Sec. 5 we describe results obtained from synthesis; in Sec. 6 we investigate the effects of 
unresolved magnetic features; in Sec. 7 we discuss the results and draw our conclusions. 



2. Data 

We investigated the behavior of the two Fe I lines at 617.3 and 630.2 nm in data 
acquired at the NSO Dunn Solar Telescope with IBIS. In order to interpret the results 
derived from observations in terms of properties of the observed solar plasma, we also 
synthetized the two lines in semi-empirical atmospheric models representative of quiet Sun 
and magnetic regions. 



2.1. Observations and data analysis 

The IBIS data set employed for this study was taken on March 17th 2011 from 14:19 
UT to 15:19 UT. It consists of 60 spectro-polarimetric scans of the two lines, each sampled 
at 22 spectral positions; the nearby O2 630.28 nm Telluric line was also sampled at 6 spectral 
positions and served for line position calibration purposes. A full scan of the two lines 
took ^ 60 seconds, so that each line was scanned in ^ 30 seconds. This temporal cadence 
is short with respect to typical evolution time scales of solar photospheric phenomena, so 
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Fig. 1. — Examples of analyzed IBIS data. Top: 617.3 nm (left) and 630.2 nm (right) 
nearby continua. Bottom: corresponding magnetic flux density maps derived with the COG 
method. The FOV is approximately 40 x 90 arcsec. 
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that the integrity of each hne during a scan can be assumed to be preserved. The Field of 
View (FOV) was approximately 45x90 arcsec and included a facular region and a small 
pore located at [16.5N,10.3E] in between the leading and following parts of NOAA 11172. 
White Light (621.3 nm) frames imaging the same FOV were acquired contemporaneously to 
spectro-polarimetric data. The spatial sampling of IBIS and WL frames was 0.17"pixel^^. 



20041 ) was tracking on 



The High Order Adaptive Optics System of the DST (iRimmele 
a small pore in the FOV during the acquisition, thus improving image quality and stability. 
The data set was reduced and compensated for instru mental blue-shift, and for instrumen t 
t elescope polarizati on by standard procedures (i.e 



and 



2009 



Cauzzi et al 



2008 



Viticchie et al 



Judge et al. 



2010l ). Residual aberrations introduced by the at mospheric seeing 



were 



20061) 



estimated by applying Multi Frame Blind deconvolution technique ( iVan Noort et al 
to WL images and were employed to destretch the spectro-polarimetric images. The 
residual crosstalk between the Stokes parameters and the rms noise level were estimated to 
be 2 • 10~^ and of the order of 10^^ the intensity in the nearby continuum, respectively. An 
infrared blocker in the optical path re duced the spectral str ay-light contamination, so that 



its residual, estimated as described in 



Criscuoli et al. 



(120 llh is below 8 % and 4 % for 617.3 



nm and 630.2 nm, respectively. For the purpose of this study, we analyzed a subsample of 
10 scans having the highest rms-contrast in WL images. Figure [U shows an example of the 
analysed data. 

We estimated the magnetic flux de nsity along the Line- of-Sight in the FOV by applying 



the Center Of Gravity method ( COG 



left- circularly polarized signals. 



Rees fc Semel 



19791 ) to the measured right- and 



Uitenbroekl (120031 ) showed that the COG provides good 



estimates of the magnetic flux at the formation height of typical photospheric line core 
for fields lower than 1.5 kCauss. To estimate uncertainties of magnetic flux density values 
derived with the COG we applied standard error propagation to the COG formula. We 



- 8 - 



found that the uncertainty is on average equal to 8 Gauss for flux density values smaller 
than 150 Gauss; it increases monotonically for larger flux density values, although the 
relative uncertainty is approximately constant and equal to 3%. The scatter of values 
increases with the decrease of the estimated flux. The average value of the uncertainties 
computed over the whole FOV is 9 Gauss, with a standard deviation of 6 Gauss. Note, 
however, that the finite polarimetric sensitivity of IBIS sets a lower limit of confidence on 
magnetic flux density estimates. To estimate this limit we computed the average absolute 
flux over those pixels in the FOV whose total polarization signal [l 

Next, in order to characterize the shapes of the lines, we computed the core intensity 
(IC, hereafter), the equivalent width (EQW, hereafter), and the full width at half maximum 
(FWHM, hereafter) of the two lines at each image pixel of the FOV. In particular, the IC 
and FWHM were computed by fitting the line sampling with a Gaussian function. We 
computed the three line shape parameters only for those pixels in the FOV with a magnetic 
flux below 1000 Gauss,thus excluding the small pore and the area around it from the 
analysis, but including network and facular elements. This allowed us to investigate the 
dependence of line shape parameters on magnetic flux density variations. 



^The total polarization signal is defined as / \Ju{Xf + V{Xf + Q{X)^dX/ Ic, where Ic is 
the continuum intensty of Stokes I, and U, V and Q are the corresponding Stokes parameter, 
is lower than 3 mA and found a value of 20 Gauss with a standard deviation of 10 Gauss. 
We therefore defined as quiet- Sun all those pixels whose absolute magnetic flux density is 
smaller that 30 Gauss (that is, below one standard deviation from the mean). 
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2.2. Spectral synthesis 



We performed synthesis of the two Fe I line profiles in various one- dimensional 
semi-empirical static models representative of various thermodynamic conditions of the 



plasma and for different values of m agnetic field strengt 



FAL-C, FAL-E, FAL-H and FAL-P flFontenla et al. 



IS. In particular, we employed 



19991 ) models, representative of quiet 



Sun, network, facular and bright facular regions, respectively. The spectral syri t 



nesis 



was performed in local thermal equilibrium with the RH code (lUitenbroek 



2002 . 



20031 ). 



The Lande factor is 2.5 for both lines; the other atomic data adopted for the synthesis 
are summarized in Table [1] The three line parameters IC, FWHM and EQW were then 
computed for emergent spectra calculated at vertical line of sight. 



3. Response Functions 

In order to discuss the results obtained from the two Fe I lines we first 
investigated their formation heights. To this aim we analysed results derived 
from spectral synthesis of the two lines through FAL-C and FAL-P models, the 
latter with uniform magnetic field of 1000 Gauss. We then computed the Response 
Function (RF, hereafter) of IC of both lines to perturbations of either the temperature 
(T_RF) or magnetic field (B_RF) on both atmosphere models. 

Figure [2] shows the T_RF of IC for both lines and models. We find that the 617.3 

Line (nm) loggf/ x (sV) 
617.3 -2.88 2.22 
630.2 -1.24 3.68 



Table 1: Lines Atomic Data 
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Fig. 2. — Temperature Response Functions of line core intensity obtained in FAL-C (contin- 
uous) and FAL-P (dashed) for the 630.2 nm (black) and the 617.3 nm (red) lines. Response 
Functions are normalized to their maxima. Optical depth unity at 500 nm sets the height 
reference value. 



- 11 - 



nm line core in FAL-C forms at a wide range of photospheric heights, from to 600 
km, with a peak a t 200-300 km above optical depth unity at 500 nm (in agreement with 



Norton et al 



20061 ). The core of the 630.2 nm line instead forms in a slightly narrower 
range of photopsheric heights, from 180 to 600 km, and peaks slightly higher in the 
atmosphere at approximately 300 km. Figure [2] also shows that the core of the two lines 
forms in the deepest layers of the photosphere in magnetic features (in agreement with 



Khomenko &: CoUados 



20071 ). This is the result of the lower temperature of FAL-P with 



which is also found in computations of magnetic flux tube models (e.g^ 



Steiner 





Pizzo et al. 


1993; 


Carlsson et al. 


2004) 



20051 ) and full 3-D Magneto Hydro Dynamic simulations (e.g. 
Both computed RF^ peak at 100 km, with contribution from higher photospheric layers 
for the 630.2 nm line and from layers below optical depth unity for the 617.3 nm line. 
The bumps of both T_RF at approximately 200 km reflect the change of the temperature 
gradient with height in FAL-P model. 

Figure [3] shows the B_RF of IC for both lines in the FAL-P model. The two B_RF are 
quite similar, as both peak at approximately 190 km, with the B_RF of the 617.3 nm line 
having again a slightly larger contribution from the deepest layers than the one calculated 
for the 630.2 nm line. 

The results derived from the synthesis indicate therefore that both lines form in similar 
regions of the solar atmosphere, with the 630.2 nm line sampling slightly higher atmospheric 
heights than the 617.3 nm line. 
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4. Results from observations 

In this section we investigate the dependence of the three line shape parameters on 
magnetic flux density. 

Note that we found that the flux density values estimated from the 630.2 nm data 
are on average lower than those obtained from the 617.3 nm set, and that the discrepancy 
increases with the magnetic flux density, with a maximum of approximately 40% at 
1000 Gauss. Because, as shown in §3, the two lines form in similar regions of the solar 
atmosphere, this discrepancy is unlikely to reflect a stratification of the magnetic field with 
height. To verify this hypothesis we applied the COG to synthetic spectra of the two lines 
computed in the FAL-P model with a stratified magnetic field and obtained magnetic field 
estimates in agreement within 4%. We conclude that the discrepancy that we found in our 
data is most likely due to the Telluric blending affecting the 630.2 nm line. In the following, 
we therefore refer to fiux density estimates derived from the 617.3 nm data only. 

Results from our measurements are presented in Fig. HI Dot symbols show the 
dependence of IC, EQW, and FWHM on the magnetic fiux density for both lines, while 
asterisk symbols indicate the average value of the same line parameters computed over 
bins of magnetic fiux. Note that the bins were increasingly larger in order to guarantee 
sufiicient statistics. We found that the IC and FWHM increase with the magnetic fiux for 
both lines, whereas the EQW shows the opposite trend. All three line shape parameters 
show saturation. For both IC and EQW this saturation occurs at approximately 300 Gauss 
for both lines, while the FWHM saturates at fiux density values above 700 Gauss and 500 
Gauss for the 617.3 nm and the 630.2 nm line, respectively. Note that the maximum relative 
variations with magnetic fiux density (i.e. relative difference between values obtained at the 
largest flux and those obtained for quiet regions) are approximately 40% for IC for both 
lines, 55% and 36% for FWHM in the case of 617.3 nm and 630.2 nm respectively, and 
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-15% and -8% for EQW, respectively. 



The measurements present large scatter, especially at the lowest flux density values, 
even though the standard deviations of the values distributions in each bin (error bars in 
the plots) are rather small. The observed scatter of values reflects the sensitivity of line 
shape parameters to temperature and velocity variations of the plasma. The larger the 
magnetic flux, the more t he temperature and velocity fluctuations of convective plasma 



motions are inhibited (e.g. 



Ishikawa et al. 



2007 



Romano et al. 



20121 ) . so that the dispersion 



of the line shape parameters becomes smaller and smaller with the increase of the magnetic 
flux density. Note that the dispersion for the 617.3 nm line is larger than that obtained 
for the 630.2 nm data, especially when comparing IC. This is due to the lower excitation 
potential of the 617.3 nm line, which makes the co re intensity of this line more sensitive to 



temperature and velocity fluctuations (iGray 



1992 



, p. 285-287). 



The IC in particular depends on both temperature and magnetic field in a manner 
that is difficult to disentangle. Results from numerical models and simulations of magnetic 
flux tubes show in fact that the photospheric temperature within a tube rises following an 
increase of magnetic field , due to the inten sification of the lateral heating effects from the 



surrounding plasma (e.g. 



Pizzo et al. 



19931 ). As also qualitatively shown by results 
obtained from numerical synthesis reported in Sec. [5l an increase in either of these 
quantities, the field strength, via the Zeeman effect, or the temperature, via irradiation, 
leads to an increase of IC, so that it is difficult to discern between the two effects. 

The trend with the magnetic flux density of FWHM found for both lines suggests 
instead that this line shape parameter is more sensitive to magnetic field variations and 
velocity effects than to temperature. In fact, from the rise of temperature within flux tubes 
mentioned above, we would expect the FWHM to decrease with the magnetic flux. This 
decrease is compensated on one hand by the fact that magnetic flux concentrations tend 
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to be located in or nearby downflow regions, which are characterized by shallower and 
broader profiles, and on the other hand, to the increase of the splitting of the lines with 
the magnetic flux. This latter effect becomes more and more important with the increase 
of the field (and therefore of the flux), due to the inhibition of convective flows mentioned 
above. This interpretation is confirmed by results obtained from numerical synthesis of line 
profiles discussed in Sec. Ei 

We also note that at low magnetic flux density values, the FWHM of 617.3 nm line 
increases slower with the magnetic flux density than for the 630.2 nm. This indicates that 
for the former line the effect of temperature counteracting Zeeman broadening is larger than 
for the 630.2 nm. Note that the maximum difference between the FWHM at saturation 
and the FWHM of average quiet Sun agree within the errors for the two lines (i.e. is 
approximately 0.06 A), as expected from Zeeman splitting theory for lines with equal Lande 
factor. 

Finally, the dependence of EQW on the magnetic flux reflects variations of both the 
width and depth of the studied lines. The rapid saturation of EQW with the magnetic 
flux is due to the fact that the line core intensity increases more rapidly than the FWHM 
for a given increase of magnetic flux. As formally shown in the Appendix [Al 
and confirmed by results obtained from numerical synthesis shown in Sec. [S], 
the EQW of the two iron lines is almost independent of the magnetic field and is mostly 
sensitive to temperature variations. Figure H] therefore suggests that the magnetic features 
in our data are characterized by similar temperatures. 
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5. Results from numerical synthesis 

In previous section we liave analysed tlie dependence of tlie tliree shape parameters of 
the two iron lines on the magnetic flux density. As already noticed, this dependence is not 
explicit, because of the indirect effects of the magnetic fleld on other physical parameters. 
In particular, the indirect effect on the temperature stratification highly affects the shape 
of a line. In order to separately investigate the effects of magnetic field and temperature on 
line shape parameters we have therefore analysed results obtained from numerical synthesis 
through various one dimensional models, as described in Sec. 12.21 Results are presented in 
this section. 

To investigate the effects of temperature on the shapes of line profiles we have 
compared results obtained from the various FAL models, which are characterized by 
different temperature stratifications. Plots in Fig. [5] show the relative differences of the line 
parameters computed for FAL-E, FAL-H and FAL-P models with respect to line parameter 
obtained with FAL-C model. Since the three models are characterized by increasing 
temperature at the formation heights of the two lines, the plots show that the increase of 
temperature causes IC to increase and the FWHM and EQW to decrease. While for IC and 
EQW we found changes up to 50% and more, the FWHM changes do not exceed 10%. 

To investigate the effects of the magnetic field, we considered line profiles obtained with 
FAL-C model in the case of fields constant with geometric height and whose strength varied 
from 50 to 1000 Gauss. Results reported in Fig. |6] show the relative difference between line 
parameters obtained for the different field strength values with respect to those obtained in 
the case of no magnetic field. We note that the increase of the magnetic field strength leads 
to a large increase of IC and FWHM (more than 100% in the case of IC), especially at the 
largest magnetic strength values, while EQW stays fairly constant (variations are smaller 
than 1% for both lines), thus confirming the analytical derivation presented in Appendix Rl 
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Both plots show that the two hnes are differently affected by changes of plasma physical 
properties, mainly because of their different atomic parameters. The plot in Fig. [2] shows 
for instance that the difference of the shapes of the T_RFs obtained for FAL-C and FAL-P 
models is larger for the 617.3 nm line than for the 630.2 nm line. In particular, while the 
T_RFs of the 630.2 nm line obtained with the two models appear mostly shifted in depth 
respect to each other, the two T_RFs obtained for the 617.3 nm line also show different 
shapes. The 617.3 nm line therefore experiences larger changes due to temperature effects 
than the 630.2 nm. 

A direct comparison of plots shown in Fig. O and in Fig. [H] is not straightforward, 
as they show the variations of the line parameters with respect to two different physical 
quantities, the temperature and the magnetic field strength, respectively. Nevertheless, 
results from magnetic flux tube models show that within magnetic flux tubes in radiative 
equilibrium the temperature at the formation heigths of the two investigated lines is 
positively correlated with the field strength, so that, for instance, model FAL-P can 
reasonably have associated a field strength of 1000-1500 Gauss. This allows a comparison of 
the trends of the line parameters with respect to the two physical quantities. The plots show 
that IC increases with the increase of both temperature and field strenght, in a manner 
difficult to disentangle. Opposite trends are instead found in the case of the FWHM, as 
an increase of temperature makes this line parameter to decrease while an increase of 
magnetic field strength makes it to increase. Comparison of the value of FWHM variation 
obtained for FAL-P with the value obtained for 1000 Gauss show that magnetic field effects 
are more likely to dominate. Variations of the EQW found in the case of the increase of 
the temperature are more than one order of magnitude larger than those obtained when 
increasing the magnetic field strength, again confirming that the direct effect of Zeeman 
splitting on this line parameter is negligible for the two lines. 
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Finally, it is important to notice that variations found from the numerical synthesis 
presented here are meant to quahtatively interpret observational results presented in 
Sec. m as they are representative of the ideal case of unitary magnetic filling factor and 
observational effects like finite spatial and spectral resolution and scattered light have not 
been taken into account. 



6. Effects of unresolved magnetic features on observed lines 

In this section we investigate whether and to what extent, line shape parameters are 
affected by small amounts of magnetic flux not resolved in the observations. 

To this purpose we analysed average line shape parameters obtained in images that 
were synthetized by replacing line profiles of pixels whose magnetic flux density exceeds 
a certain threshold with the average quiet-Sun profile. By varying the magnetic flux 
threshold, images of successively lower average magnetic flux density (in particular from 
approximately 90 to 10 Gauss) were constructed. Note that the quiet Sun profile was 
computed averaging profiles of those pixels of the original images whose measured magnetic 
flux density is lower than 30 Gauss (see Sec. 2). The line shape parameters were then 
determined from the average profiles of each of these synthetic maps, thus mimicing 
observations that do not resolve the remaining magnetic features. Figure [7] shows the 
relative differences between line shape properties from the synthetic frames (characterized 
by different values of the average magnetic flux) and those of the synthetic frame with the 
lowest magnetic flux density (i.e. 10 Gauss). In the case of IC and FWHM the error bars 
are derived from the uncertainties of the fitting of the l ine to a Gaus sian function (see Sec. 



Cayrell (1l988l ) formula. We found 



2), while for the EQW the errors are estimated by the 
that the average line shape parameters vary by a few percents with the increase of the 
average flux density. The maximum variation that we found are approximately 6% in the 
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case of IC for both lines, 4% and 2.5% in the case of FWHM for the 630.2 nm and the 
617.3 nm, respectively, and -1.8% and -3.5% in the case of EQW. Variations are larger in 
617.3 nm in the case of EQW, as a result of the already mentioned larger sensitivity of this 
line to temperature. For the same reason, larger variations of FWHM with the magnetic 
flux density are instead found for the 630.2 nm. As explained above, in fact, for the 617.3 
nm line the effect of temperature counteracting the Zeeman line broadening is larger, thus 
reducing the increase of the FWHM of this line with the flux. 



7. Discussion and Conclusions 

The measured line shape parameters of the two Fe I lines vary with magnetic flux 
because of the direct effect of the Zeeman splitting and the indirect effect of temperature 
variations induced by suppression of convective motions. IC and FWHM increase with 
the magnetic flux, while the opposite trend was found for EQW (Fig. H]). Results from 
numerical synthesis show that in the case of IC, the trends reflect both the Zeeman splitting 
and the temperature effects in a manner that is difficult to disentangle, whereas the FWHM 
is mostly sensitive to the Zeeman effect, in particular for larger field strenght. We have 
also shown that the dependence on magnetic field of the EQW is negligible at the values 
of magnetic flux investigated, so that the observed trend must be ascribed to temperature 
variations induced by the presence of magnetic field concentrations. 

We also found that the three line shape parameters show saturation with the magnetic 
flux density (Fig. H]). Interpretation of these results in terms of magnetic flux tube models 
suggests that features characterized by small magnetic flux values are mostly spatially 
unresolved in our observations and that the pixels characterized by the largest magnetic 
flux values (above the saturation threshold) belong to spatially resolved features that may 
result from the clustering of magnetic elements. According to flux tube models, in fact, an 
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increase of the magnetic field leads to an increase of the temperature within the tube in 
line-forming layers. Instead of a saturation we should therefore observe IC and FWHM to 
monotonically increase and EQW to monotonically decrease with the flux. In particular, 
the saturation of EQW suggests that the magnetic features that we observed are mostly 



characterized by very similar photos pheric temperatu res 



clusters of magnetic flux tubes (e.g. 



Criscuoli &: Rast 



and is in agreer nent with results recently rep orted by 



notice that also 



Harvey fc Livingston! ( Il969a 



2009 



lis is expected from models o f 



Fabiani Bendicho et al. 



Viticchie et al 



1992) 



mm . It is worth to 



b|) concluded that the proflles of the Fe I 525.0 



and 525.3 nm lines suggest the "quantization" of the magnetic flux in the photosphere. 

We also investigated the effects of unresolved magnetic fleld on the shapes of line 
profiles. Results reported in Fig. [7] show that variations of a few percent in the average 
shapes of the investigated lines are produced for modest (below 100 Gauss) values of average 
magnetic flux. Note that the average flux values investigated are compatible with 
those err iployed in irradiance reconstructio n techniques to identify quiet Sun 



regions ( jFligge et al.l 



2000l : iBall et al. 



2011 



, and references therein). In the case of 



FWHM, larger variations with the average magnetic flux are found for the 630.2 nm than 
for the 617.3 nm, whereas in the case of EQW larger variations are found for the 617.3 nm. 
The 630.2 nm line is therefore more favorable to investigate magnetic fleld variation, while 
the 617.3 nm is more suitable to investigate effects of temperature variations. 

It is interesting to notice that our estim ates of EQW variations are in agreement 
with values reported by 



Fabbian et al. 



mm for the Fe I 608.27 nm and 624.07 nm lines. 



These authors analyzed snapshots from Magneto Hydrodynamic Simulations to investigate 
the variations of average EQW values of some iron lines with the increase of the average 
magnetic flux in the simulations. Comparison of our Fig. [7] with their Fig. 4 reveals that 
the agreement is better for the 617.3 nm line, as this has an excitation potential similar 
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to that of 608.27 nm and 624.07 nm lines. Their plots also indicate that the variations of 
EQW mostly reflect variations of the average temperature in their simulations. The small 
dependence in their plots on the average magnetic flux is due to the wavelenght shift of 
the opacity discussed in the Appendix. For completeness, we notice that those authors 
also reported a large dependence of EQW of the Fe I 1.5 yum line on the magnetic flux. 
This is due on one hand to the fact that this line forms in regions where the variation of 
temperature induced by the magnetic field is negligible, on the other to its large Lande 
factor so that this line is not in the weak field regime. 

The amount of variations of line shape parameters due to the presence of unresolved 
magnetic flux that we found is measurable with instruments like HINODE/SOT, 
SOLIS/VSM and SDO/HMI. Line shape variations of the investigated lines observed at disk 
center can be employed for studies of long term temporal variations of physical properties 
of quiet Sun. In particular, the results presented here show that it would be possible to 
disentangle thermal and magnetic effects by comparing variations of different properties of 
the two Fe I photospheric lines. For instance, a decrease over the magnetic cycle of EQW 
together with an increase of FWHM would indicate variations of quiet Sun temperature 
induced by the presence of unresolved magnetic features. On the other hand, a decrease of 
the EQW together with a decrease of the FWHM would indicate an increase of the effective 
temperature of the plasma. 



Finally, it is worth to notice that 



Martinez Gonzalez 



showed that the same 



profiles of Fe I 630.1 and 630.2 nm lines can be produced by different thermodynamic 
stratifications, magnetic field and microturbulence values, so that results obtained from 
spectro-polarimetric inversions should be treated with caution. However, the method that 
we propose is less affected by this degeneracy. In fact, we have shown that the EQW 
parameter of the two iron lines is almost insensitive to magnetic field effects, while the 



FWHM is mostly sensitive to magnetic field effects. The EQW variations can therefore 
be utilized to derive information on temperature variations, while the FWHM can be 
utilized to derive information on magnetic flux variations. No te, however, that a proper 



modeling of the observed data (e.g. 



Holzreuter fc Solanki 



2OI2I ) is required to quantitatively 



estimate these variations. In particular, modifications of the velocity distributions, which 
can result from either variations of the magnetic flux or of the effective temperature, can 
also affect the shapes of lines. Nevertheless, since velocity distributions are known to 



narro w in the case of increase of magnetic flux (e.g. 



Ishikawa et al. 



2007 



Morinaga et al 



20081 ). while they broaden and skew toward upfio^y velocities in the case of an increase of 



the effective temperature (e.g. 



Gray fc Toner 



1992 : 



Nordlund fc Dravins 



19901 ) ■ we expect 



velocity effects not to hamper the possbility of disentangling magnetic flux and effective 
temperature variations from observations of properties of the two lines. We plan to 
investigate these effects in the future by analysing results from numerical simulations of the 
solar photosphere. 
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A. Dependence of Equivalent Width on magnetic field: analytical derivation 

Here we investigate the dependence of the EQW of a magnetically sensitive spectral 
line on the magnetic field. It is common practice to express the emerg ent line profile /(A) 



in the presence of Zeeman splitting as the sum of three Gaussians (e.g. 



Borrero 



20081 ) : 



,,,, _ (l-/o)sin^7 (^-^o)^ (1 - Jo)(cos^7 + 1) , - „ , ^- - . , , 
J(A) = 1-^ 1 le ^ ^^^^ ^ ^x(e ^ — +e ^ — ) (Al) 

The two Gaussian functions shifted by the amount Aq ± As represent the left- and right-hand 
circularly polarized a profiles, while the Gaussian function centered at Aq represents the 
profile of the linearly polarized vr component. Jq and AA are the core intensity and Doppler 
width of the line, in the absence of magnetic field; 7 is the inclination angle with respect to 
the line of sight; \b = CgejfX^B , where C is a constant, Qejf is the effective Lande factor, 
Ao is the central wavelength of the line, and B is the magnetic field strength. 

In the case of magnetic field oriented along the line of sight (7 = 0) the previous 
formula simplifies to: 

. (1 - -^0) / (A-Ao-Afl)2 (A-Ao+Afl)2 

/(A) = l-^ — —^x{e Ka2 +e Xa2 ) (A2) 
The EQW of the line in that case is, by definition: 



EQW:=^-^[ re-'-''^dX+ H e'^'^^dX] (A3) 
/o Jo 



2 



The two terms in the sum are equal as they differ only for by a shift within the infinite 
domain. Since the two terms represent the EQWs of the left {EQW~) and right {EQW~^) 
circular polarized signal, then EQW = EQW" = EQW^ (note that it is easy to verify 
that these equalities hold for any shape of the profiles, as far as the two circular polarized 
signals are symmetric). 

We want now to study the variation of the EQW with the magnetic field. To this aim 
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we differentiate Eq JA3l with respect to B: 



dEQW 
dB 



In d 



2 dB 



e d\ + 



(A-A(i + Ag)^ 

e ^3:^ (iA 



(A4) 



Making the variable change A = A — Aq, and changing the order of integration and 
differentiation, we obtain: 



dEQW 
dB 



[I - Io)Cgeff\l 



-^dX + 



(A + Ab) _(^±^,r- 

-^dX 



(A5) 



where we have assumed that the rest wavelength is much larger then the Zeeman splitting. 
It is easy to see that the right-hand side of EQ. IA5I is null, as the terms in the integrand 
that are proportional to A are uneven, and the terms that are proportional to As cancel. 
Thus, under the assumption that Iq and AA do not depend on the field, the EQW of a 
Zeeman split line is independent of the field strength. Howevere, this assumption does not 
hold in a stratified atmosphere, because the Zeeman splitting changes the opacity at a 
given wavelength, and therefore alters the height and physical p arameters the line samples. 



Indeed, as shown by previous studies (e.g. 



Fabbian et al. 



20101), the EQW of a magnetic 



sensitive line is sensitive to the strength of the magnetic field. These variations are small 
in the weak field case (i.e. Aq ± A^ < AA), and are non-negligible only for large fields. 
Moreover, the sensitivity of the EQW to the field depends on the atomic parameters of the 
line, as these determine its opacity. 
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Fig. 3. — Magnetic Field Response Functions of line core intensity obtained in FAL-P, 
representing facular regions, for the 630.2 nm (black) and the 617.3 nm (red) lines. 
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Fig. 4. — Variations of line shapes parameters IC, FWHM and EQW for 617.3 nm (top) and 
630.2 nm (bottom) lines with the magnetic flux density. 
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Fig. 5. — Relative variation of line parameters obtained from spectral synthesis for different 
atmosphere models with respect to FAL-C. Diamond: 617.3 nm; Triangles: 630.2 nm. 
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Fig. 6. — Relative variation of line parameters with the magnetic field strength obtained 
from spectral synthesis in FAL-C assuming uniform magnetic field stratification. Diamond: 
617.3 nm; Triangles: 630.2 nm. 
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Fig. 7. — Relative variation of average line shape parameters with the average magnetic flux 
density computed on synthetic frames (see text). Diamond: 617.3 nm; Triangles: 630.2 nm. 



